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Abstract. In the framework of the magnetar model for the Soft Gamma 
Repeaters and Anomalous X-ray Pulsars, we consider the structure of 
neutron star magnetospheres threaded by large-scale electrical currents. 
We construct self-similar, force- free equilibria under the assumption of 
axisymmetry and a power law dependence of magnetic field on radius, 
B oc r~'^'^~^P' . A large-scale twist of the field lines softens the radial de- 
pendence to p < 1, thereby accelerating the spindown torque with respect 
to a vacuum dipole. A magnetosphere with a strong twist {B^/Bq = 0(1) 
at the equator) has an optical depth ~ 1 to resonant cyclotron scattering, 
independent of frequency (radius), surface magnetic field strength, or the 
charge/mass ratio of the scattering charge. We investigate the effects of 
the resonant Compton scattering by the charge carriers (both electrons 
and ions) on the emergent X-ray spectra and pulse profiles. 



1. Nature of the AXPs and SGRs 

The Anomalous X-ray Pulsars and Soft Gamma Repeaters are neutron stars 

which share similar spin periods, P = 5 12 s, characteristic ages, P/P = 

3 X 10'^ — 4 X 10^ yr, and X-ray luminosities, Lx = 3 x 10^^ — 10^^ erg s~^, well in 
excess of the spin-down luminosities (Thompson et al. 2001). They are almost 
certainly young and isolated: no evidence for a binary stellar companion has yet 
been detected in any of these sources, and a few are convincingly associated with 
young supernova remnants. The overlap between the SGR and AXP sources in 
a three-dimensional parameter space (P, P and Lx), and the observed variabil- 
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ity in their X-ray output, provides circumstantial evidence that they share a 
common energy source: the decay of a very strong 10^^ G) magnetic field. 

The persistent emission of the SGRs and AXPs has both a thermal and a 
non-thermal component. In the case of the SGRs, this emission becomes brighter 
after periods of bursting activity and involves a comparable release of energy to 
the outbursts (averaged over time). 

2. Twisted Neutron Star Magnetospheres 

The magnetic fields of magnetars are most likely generated by a hydromagnetic 
dynamo as the star is born, and may be associated with rapid initial rotation 
(Duncan &: Thompson 1992). A strong twist in a ~ 10^^ G magnetic field will 
relax at intervals as the field is transported through the deep interior of the 
neutron star. Even if the electrical current were initially confined to interior 
of the star, the Lorentz force would become strong enough to deform its crust, 
thereby twisting up the external magnetic field. In such a situation, a persistent 
current will flow through the magnetosphere, supported by emission of light ions 
(e.g. H, He, C) and electrons from the neutron star surface. The decay of this 
current outside of the star is an efficient mechanism for converting magnetic 
energy to X-rays, and for inducing rapid variations in the X-ray flux. 

To see how the persistent current will modify the structure of the mag- 
netosphere, we find axisymmetric force-free equilibria outside a (non-rotating) 
spherical surface, V x B = a{V)'B. These equilibria form a one-dimensional 
sequence labeled by the flux parameter V = V{R,9), with poloidal magnetic 
field Bp = W X ^/RsinO. As a major simplification we consider self-similar 
configurations 

V = VoiR/Rmr^Fie), BiR, 9) ~ F{9) x (iJ/i^Ns)-^'^^^ (1) 

following Lynden-Bell & Boily (1994). The radial index p is uniquely determined 
by a single parameter C, which is related to the strength of the current: 

p{p + 1)F + sin2 9 J^^ = -CF^+yP. (2) 
a(cos 9y 

Choosing a dipole field at a zero current, p{C) is determined by the three bound- 
ary conditions, Br{R,9 = 7r/2) = 0, Br{RNS,9 = 0) = const, B^{R,9 = 0) = 0. 
The index p is most conviently expressed as a function of the net twist A^n-S 
between the north and south magnetic poles (Fig. la). 



3. Resonant Scattering 

The current-currying charges also provide a significant optical depth to resonant 
cyclotron scattering. For a particle of charge Ze and mass M, the resonant cross- 
section is aj-esiuj) = {ir'^Ze^ /Mc) (1 -|- cos^ 9kB)^{^ — ^c)- The optical depth is 
determined by relating the particle density nz to the twist in the magnetic field 
through {Ze)nzvz = ez{c/4Tr)\V x B|. In our self-similar model, 

("^) '^'■''^ ^ ("^ ^^^^ ^^^^ 
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Figure 1. (a) The radial index p as a function of the twist A0n-s- 
(b) The actual polar field inferred from spin parameters P and P, 
compared with the magnetic dipole formula. 

Thus the optical depth is Tj-cs ^ {"Vz/c)^^ for strong twists A</)n-s ~ 1- Re- 
markably it is independent of the mass and charge of the scatterers, the radius, 
or the resonant frequency. 



4. Implications for X-ray Spectra and Pulse Profiles 

• Surface Heating. If ions and electrons supply the current, then the stellar 
surface is heated at the rate 
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(The electrons are electrostatically accelerated downward above the an- 
ode.) This is comparable to the observed luminosities of AXPs and SGRs. 
A global twist will decay on the timescale 
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Resonant Comptonization. Both ions and electrons can be expected to 
move mildly rclativistically where they resonantly scatter 1 — 10 kcV pho- 
tons. The product of the mean frequency shift {Alo/lo) per scattering 
with the expected number of scatterings is 0(1) when A(?!)n-s ~ 1 radian. 
Therefore, multiple scattering at the cyclotron resonance by the moving 
charge carriers will create a non-thermal spectral tail to the X-ray flux 
emerging from the surface. At optical depths <^ 1, the hardness of the 
spectrum increases with the number of scattering and thus with the reso- 
nant optical depth. 

Pulse Profiles. The emergent pulse profile is strongly modified by magne- 
tospheric scattering. Three effects enter here: the optical depth Tres{d) 
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is anisotropic, vanishing toward the magnetic poles; the resonant sur- 
face is aphserical; and the scattered radiation tends to be beamed along 
the magnetic field (in part because of the motion of the charge carriers). 
The pulse profile will be approximately frequency-independent in this self- 
similar model. 

• Ion Cyclotron Resonance. The ion component of the current will generate 
a comparable optical depth to the electron component in a self-similar mag- 
netosphere (at frequencies below the surface cyclotron frequency). Since 
the ion cyclotron resonance sits much closer to the star (at ~ 10—20 km for 
2 — 10 keV photons), it is more sensitive to the presence of higher magnetic 
multiples. Thus, a combination of ion and electron scattering in a multi- 
polar magnetic field will produce an energy dependent total profile which 
is more complex at higher energies. An ion cyclotron emission line would 
be a clear observational signature of surface heating by magnetospheric 
charges. 

• SGR/AXP spindown. When twisted, the magnetic field drops off more 
slowly than ~ R~^. Thus the real polar field inferred from P, P is smaller 
than the magnetic dipole formula would imply, while the braking index 
becomes n = 2p + 1 < 3 (Fig. lb). Our model predicts that the spindown 
rate increases with the optical depth to resonant scattering, and hence with 
the hardness of the persistent X-ray spectrum. In fact the active SGRs 
1806-20 and 1900-1-14 both have higher P and harder X-ray spectra than 
any AXP. The quiescent SGR 0525-66 has a softer spectrum. We predict 
that its spindown rate, when measured, will be intermediate between these 
sources and the AXPs. 

• Giant Flare Mechanism There are two generic possibilities for the produc- 
tion of the giant flares of the SGRs, in the framework of our model: (i) 
giant flares may result from a sudden change (unwinding) in the internal 
magnetic held, implanting a twist into the magnetosphere; or alternatively 
(ii) the twist may build up more gradually in the magnetosphere, leading 
to a sudden relaxation in close analogy with Solar flares. Measurements 
of hardening/softening of the persistent X-ray spectrum, changes in pulse 
profile, and changes in spindown rate provide a way of discriminating be- 
tween these possibilities. In particular, the simplified pulse profile observed 
in SGR 1900+14 after the 27 August 1998 giant fiare (Woods et al. 2000) 
can be explained by an increase in the currrent at the radius of the electron 
cyclotron resonance. 
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